Symmetry, Integrability and Geometry: Methods and Applications 



Vol. 1 (2005), Paper 018, 15 pages 



in 



o 



d 



o 



> 
X 



Ermakov's Superintegrable Toy 
and Nonlocal Symmetries 

P.G.L. LEACH ^ , A. KARASU (KALKANLI)^ , M.C NUCCI? and K. ANDRIOPOULOS^ 

t School of Mathematical Sciences, Howard College, University of KwaZulu- Natal, 
Durban 4041, Republic of South Africa 
E-mail: leachp@u kzn.ac.za\ 



t Department of Physics, Middle East Technical University, 06531 Ankara, Turkey 



E-mail: akarasu@metu.edu.tr 



t Dipartimento di Mathematica e Informatica, Universita di Perugia, 06123 Perugia, Italy 
\ E-mail: \nucci@unipg. i1\ 

' Department of Information and Communication Systems Engineering, 
^ . University of the Aegean, Karlovassi 83 200, Greece 

E-mail: kand@aegean.gr 



Received September 19, 2005, in final form November 11, 2005; Published online November 15, 2005 
Original article is available at http://www.emis.de/journals/SIGMA/2005/Paper018/ 



Abstract. We investigate the symmetry properties of a pair of Ermakov equations. The 
system is superintegrable and yet possesses only three Lie point symmetries with the al- 
gebra sZ(2,R). The number of point symmetries is insufficient and the algebra unsuitable 
£>. ' for the complete specification of the system. We use the method of reduction of order 

, to reduce the nonlinear fourth-order system to a third-order system comprising a linear 

' second-order equation and a conservation law. We obtain the representation of the com- 

plete symmetry group from this system. Four of the required symmetries are nonlocal and 
the algebra is the direct sum of a one-dimensional Abclian algebra with the semidirect sum 
ly-j , of a two-dimensional solvable algebra with a two-dimensional Abelian algebra. The prob- 

ed ' lem illustrates the difficulties which can arise in very elementary systems. Our treatment 

demonstrates the existence of possible routes to overcome these problems in a systematic 
fashion. 
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1 Introduction 

In 1880 V. Ermakov j2| introduced a system of second-order ordinary differential equations, 
videlicet 

x + uj 2 (t)x = — ^, (la) 



x 



y + L0 2 {t)y = 0, (lb) 

for which an invariant could be constructed by means of the operations x • Ijlbj) — y • (|laj) . 
multiplication by an integrating factor xy — xy, and integration with respect to time. The 
invariant obtained is 



(xy - xy) 2 + ( - 
\x 



2 

X 



(2) 
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In the context of Ermakov's system the dependent variables x and y could be interpreted as 
spatial variables and the independent variable t as time although this is not necessary. In such an 
interpretation the Ermakov system is a pair of coupled oscillators, one of them being nonlinear, 
with time-dependent frequency uj 2 (t). The Ermakov invariant is then a generalisation of the 
conservation of angular momentum |23| . 

In another context and era the motion of a charged particle in an axially symmetric electro- 
magnetic field was modelled in the linear approximation as the time-dependent linear oscillator 

q + co 2 (t)q = 0. (3) 

By an application of Kruskal's asymptotic method ^§1 Lewis [251 EB HZ] obtained an invariant 



{pp-pq) + 



, 2' 

2 , / 9 



in which p(t) is a solution of the auxiliary equation 



(4) 



p + u ?(t)p = L (5) 



and p is the momentum canonically conjugate to q. Lewis and Riesenfeld [2B] provided a quan- 
tum mechanical treatment of the time-dependent linear oscillator through the invariant (jlj as 
a linear operator in the context of the Schrodinger Equation. This implied connection of I 
with an Hamiltonian operator was confirmed by Leach [2Uj who related the Hamiltonian of the 
time-dependent linear oscillator 

with I through a time-dependent linear generalised canonical transformation [HI Cj 

Q = l, P = pp-pq, T = [ p~ 2 (t)dt, 
P J 

where p(t) is again given by (JSJ). 

The approaches of Ermakov and Lewis are different even though they embrace the same three 
mathematical objects. The order of introduction is different. Ermakov has a system, which could 
in principle be coupled although this is not found in the original treatment, and the system leads 
to an invariant. Lewis has a single equation and an invariant which relies upon the existence of 
a second equation. The modes of interpretation, the Ermakov invariant as a generalisation of 
the conservation of the magnitude of angular momentum and the Lewis invariant as a type of 
energylike Hamiltonian, reflect the differences in approach 1 . 

A further complication had already been provided by Pinney |46U37) who related the solutions 
of the auxiliary equation © to that of the time-dependent linear oscillator JHJ) according to 

l 

p(t) = (Au 2 + 2Buv + Cv 2 ) 2 , (6) 

where u(t) and v (t) are any two linearly independent solutions of (JHJ and the three constants 
A, B and C are related according to AC — B 2 = 1/W 2 , where W := uv — uv is the constant 
Wronskian of the basis of solutions of © . 

1 In this respect the interpretation by Eliezer and Gray |S| of the Lewis invariant as a generalisation of angular 
momentum was rather prescient in that it preceeded by a few years the spread of knowledge of Ermakov systems 
beyond the region of Kiev. That Giinther et al 1121 could demonstrate the failure of the interpretation when the 
number of spatial dimensions was increased underlies the fortuitous nature of the interpretation. 
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More recent literature |3Ul 1331 |5] links the solution © to that of the linear third-order 
differential equation 

y +4uJ 2 (t)y + Acu(t)u{t)y = 0, (7) 

where y = ^p 2 , and its integration by means of the integrating factor y. This connection becomes 
even more manifest when one recognises that (0 is the equation satisfied by the coefficient 
functions of the three generators of the algebra sl(2,M) possessed by all linear nth-order ordinary 
differential equations of maximal symmetry 2 and that sl(2, K) is also the algebra of the Lie point 
symmetries of the Ermakov-Pinney equation Qlaj). ©. 

Around 1978 the systems of Ermakov became known in the West - there is more than a little 
folklore about the mode of dissemination - and immediately attracted a great deal of attention; 
for a modest sampling see |471 1481 E9"] , The flood of papers did diminish somewhat, but more 
recent articles |13 (I14| IT5] indicate that the mathematical and physical riches of Ermakov systems 
have not been exhausted. 

In this paper we wish to treat a two-dimensional Ermakov system in a manner analogous to 
the treatment of the 'beloved' equation 3 

x + 3xx + x 3 = (8) 

by Karasu et al JB] i n which a doubled version of (jSJ) was studied from the viewpoint of its 
symmetries and integrability in the sense of Painleve. The collapse of the Lie point symmetry 
algebra from the single (JSJ) to its dual was dramatic. In that treatment there existed a linea- 
rising transformation to a pair of third-order ordinary differential equations and the Lie point 
symmetries of the third-order equations enabled the complete symmetry group of the system of 
second-order equations to be obtained in terms of a set of five nonlocal symmetries which was 
distinguished by being Abelian. 

The complete symmetry group of a system of ordinary differential equations is the group 
of the algebra of Lie symmetries required to specify the system completely. The concept was 
introduced by Krause |17l I18j in a study of the Kepler Problem. The absence of a sufficient 
number of Lie point symmetries for the purpose caused Krause to introduce nonlocal symmetries 
by an ingenious strategy which has proved to be useful elsewhere |41j . The necessity for the 
introduction of an Ansatz for the structure of nonlocal symmetries was successfully challenged 
by Nucci who used the method of reduction of order [35] to obtain the requisite number of 
symmetries wholly through point symmetries - naturally with a coupling to a nonlocal transfor- 
mation in the case of those symmetries which were nonlocal in the original representation. This 
same method has been found to apply to a wide variety of systems |40( I42j and has revealed 
an unanticipated algebraic identity of somewhat dissimilar problems. The present theoretical 
development of the properties of complete symmetry groups can be found in 12 02 E] . 

2 There are some variations in the precise forms of the coefficients of y and y to allow for the value of n; see 
Mahomed et al I3UI . Indeed the third-order equation Q as adjusted for the order of the equation under study, 
provides a unifying feature to the study of the Lie algebras of nth-order ordinary differential equations as it is 
always present for equations of maximal symmetry. For n > 3 the other symmetries come from the solution of 
an nth-order differential equation and a first-order differential equation. The solutions of the nth-order equation 
and third-order equation are related to the solutions of a second-order equation of maximal symmetry I32| . In 
fact this relationship leads to a simple demonstration of Pinney's solution (JSJ. For n = 2 another second-order 
differential equation intrudes to provide the noncartan symmetries. This second-order differential equation is the 
adjoint of the second-order equation under study |5UI . 

3l L'equation adore' as (JHJ was dubbed, perhaps a little irreverently, by the recently and sadly late Mark Feix 
in 1996 after yet another attractive property of this remarkable equation was revealed in the work of his student, 
Claude Geronimi. 
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The system which we treat is 



1 



(9a) 



x 



,3 • 



.1" 



1 



(9b) 



y 



in which each component is the simplest form of the Ermakov-Pinney equation possible and 
no coupling between the two equations exists. The reader sceptical of the potential difference 
between a single Ermakov-Pinney equation and a pair cloned as above is referred to the system 
discussed by Karasu et al ^H] i n which these differences are displayed in detail. We keep the 
system at its simplest to avoid masking the clarity of the concept under consideration behind 
a complexity of model 4 . We note that Q is related to the more usual form of Ermakov systems 
by means of a point transformation |23| . 

We examine the system © from the point of view of its Lie point symmetries and find that 
these are insufficient to specify completely the system. A recent work of Nucci et al [33] on the 
complete symmetry group of (|9a|) does provide an easy route to the complete symmetry group 
of the system © although the ease does follow from the prior passage of the route. For a single 
Ermakov-Pinney equation that route required the application of the method of reduction of 
order. The method relies upon the existence of an ignorable coordinate. In the case of @ time 
is the obvious ignorable coordinate. The choice of a new independent coordinate is not obvious 5 . 
However, a transition to polar coordinates makes the choice traditionally obvious. In polar 
coordinates the method of reduction of order leads to the combination of a linear second-order 
ordinary differential equation and a trivial first-order differential equation, i.e. a conservation 
law, reminiscent of the reduction of the MIC-Z Problem |24| in addition to problems more 
obviously connected to the Kepler Problem. 

The reduction provides a source of Lie point symmetries which can then be translated to 
the original coordinates and equations to give a suite of symmetries from which the complete 
symmetry group of system © can be determined. Fortunately the theorems developed by An- 
driopoulos et al [3] obviate the necessity for lengthy calculations involving nonlocal symmetries 
for we can use the symmetries of the complete symmetry group of the derivate system and the 
transitive properties of complete symmetry groups to arrive at the complete symmetry group 
of @- We note that in this derivation the nonuniqueness of the representation given by Nucci 
et al is replicated due to the known property of the possession by linear second-order ordinary 
differential equations of at least three representations of their complete symmetry groups |221 [2] . 

The interesting feature of system Q is that it is rather light on Lie point symmetries and 
yet, as we see below, is a superintegrable system in the standard sense in that it possesses three 
first integrals. 

The essential points of this paper fall into two parts. Individually (|9aj) and (|9bj) are endowed 
with sufficient symmetry to ensure an easy integration. In this context we mean sufficient 
Lie point symmetries as is the standard for second-order equations and usual for higher order 
equations. The simple cojoining of the two equations produces the situation in which the number 
of Lie point symmetries is formally fewer than that required for easy integration and yet the ease 
of integration remains. We use this simple system to emphasise the problems of moving from 
a scalar equation to a system of equations. The second part of the paper deals with the solution 
of the problem posed by the lack of Lie point symmetries by comparison with the individual 
scalar second-order equations. The resolution of such problems is not necessarily by means of 

4 In a sense this already happened with complete symmetry groups which entailed the use of nonlocal symmetries 
|17II18I whereas the idea can be displayed in the context of point symmetries and the most elementary of ordinary 
differential equations 

In the case of just J9al this is obviously x. 
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a recipe. Were this so, the solution of systems of nonlinear equations would be no more difficult 
than that of scalar equations. Anyone with the slightest acquaintance with the problem of the 
solution of systems of nonlinear equations knows the foolishness of that position. What we 
do in this paper reflects the two parts. Firstly we delineate the problems encountered in the 
simple act of cojoining two integrable equations as far as the usual criteria of integrability in 
the Lie theory are concerned. Secondly we demonstrate how those problems are overcome. The 
particular resolution of these problems may be unique to the system, but the techniques are 
of a general nature and belong to the standard repertoire to be used when a system displays 
refractory qualities. The techniques used here were not designed for this system. They exist and 
may be employed generally. The success of their employment depends upon the specific nature 
of the system under consideration. 

We emphasise that the system we study © is quite artificial and is studied more to illustrate 
properties and methods than to provide novel application. Yet we find it an appropriate object 
for our attentions in this 125th anniversary of its genesis in this noble city of Kiev by one of its 
now well-known sons. 



2 The superintegrable toy 

There are three Lie point symmetries of system Q, videlicet 

ri = d t , 

T 2 = td t + \{xd x +yd y ), (10) 
T 3 = t 2 d t + t(xd x + yd y ), 

with the Lie Brackets 

[r 1 ,r 2 ]=r 1 , [r 1 ,r 3 ] = 2r 2 , [r 2 ,r 3 ]=r 3 

characteristic of the algebra sZ(2,K). The system Q may be regarded as a pair of first integrals 
consequent upon the use of the integrating factors p and a, respectively, of the system 

p=0, PP=0, PP-\p 2 = \A P=\x 2 , x 3 x = A, 

a = 0, ^ aa =0, ^ aa - \a 2 = \B, ^ a= \y 2 , ^ y 3 V = B 

with the specific choice of the constants of integration A and B being unity. The third-order 
system is of maximal symmetry with 



An 


= 9 P , 


A 2 i 


= du , 


Al2 


= td p , 


A 22 


= td a , 


A13 


- h 2 d 

— 2 1 Up, 


A 23 


- h 2 d 

— 2 l U c 


Al4 


= pdp, 


A 24 


= crd a , 


Ai5 


= crdp, 


A 25 


= pda, 



(12) 



A c i = d t , 
A c2 = td t , 
A c3 = t 2 d t + 2t(pd p + ad a ) 



and only three of these survive in the integrals whether they be treated as functions or equations 6 . 

6 System li lt has the usual ambiguity of self-similar symmetries for linear systems of maximal symmetry. The 
nonlinearity of © enforces the precise form of F2 in 1101 . In terms of the variables p and a the symmetries F2 
and F3 are in the standard form for the representation of sl(2, R) appropriate to equations of the third order. The 
transformation to x and y gives the standard form for second-order equations |3U| . but we do well to recall that 
the symmetries are from the third-order system (I1H as listed in 11211 . 
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Two first integrals of system ([""} are obviously 

1 = \ \f + ^2 ) aIld J =\ [f + ^2 



The third first integral is the Ermakov invariant 

(13) 



K = \ 



(xy-xy) +^2+ ~2 



which is constructed in the same manner as for the Ermakov invariant (J""J of system (J""}. Sys- 
tem ("""J is of the fourth order and has three autonomous constants of the motion. Hence sys- 
tem (jHJ) is superintegrable. 

The autonomous integrals can all be derived from the symmetry dt if one uses Lie's method 
(cf. |2"T"|). Invariance under T\ means that the integral is autonomous, i.e. of the form F(x, y, x, y). 
The associated Lagrange's system for dF/dt = is 

dx dy dx dy 



x y x 3 y 3 



(14) 



when are taken into account. The combination of the first and third of (|14(1 leads to I 
and of the second and fourth to J. The determination of K is a little more complicated. The 
combinations y ■ (fTlk ) - x • (fTIh ) + y ■ (jTlb ) - x ■ (|TIH ) and y ■ (fTlb) - x ■ (fTlb) give 

d (xy — xy) ydx — xdy 

jE. ^_ yx — xy 

y3 x 3 

from which the third integral K follows. 

Although T\ generates the three autonomous first integrals of system Q, each integral pos- 
sesses an infinite number of Lie point symmetries. In the case of K these are [101 II 1 J 

= f(t)d t + i/(t) (xd x + yd y ) , 

where f(t) is an arbitrary, differentiable, function. In the cases of / and J the infinite classes of 
Lie point symmetries are respectively 

= d t + f(t,x,y)d y , 
Sj = d t + f(t,x,y)d x 

which are both perhaps to be expected and not as interesting as in the case of K. 

As an amusing endpoint to this discussion of the integrals we note that the integrals associated 
with T3 are 



h = hi (t± - xf + ( - 1 \, j 3 




and K as given in l)13|) • If one uses the point transformation 

T = ~, X = -, Y= y -, (15) 
t t t y J 

the integrals I and J are recovered. (Obviously K is invariant under the point transformation.) 
The transformation (|15[) is a simple instance of the generalised canonical transformation |S1 [7] 

T= [ p- 2 (t)dt, Q=l, P = PP -pq 
J P 

which has found much application in both Hamiltonian and Quantum Mechanics |2(Jl 1291 14*5] 7 . 

7 Apart from K the symmetry Y2 has the two integrals h = x(tx — x) + t/x 2 and J2 = y(ty — y) + t/y 2 . This 



Ermakov's Superintegrable Toy and Nonlocal Symmetries 



7 



3 Complete symmetry group 

The Lie algebra sl(2, M) is not found as a constituent of a complete symmetry group since 
the Lie Brackets show a redundancy of operators and this thereby leads to a violation of the 
condition of minimality which is one of the requirements for a representation of a complete 
symmetry group 144j . This then poses the problem of finding the elements of the algebra of 
the complete symmetry group for © since the only Lie point symmetries are the elements of 
sZ(2,R). Karasu et al |16j . when confronted with a similar obstacle to progress in determination, 
made use of a nonlocal transformation to a third-order system of the form of and were 
able to use the nonlocal forms of the remainder of (fT2*)l at the second-order level to express 
the complete symmetry group of the system in terms of a five-dimensional Abelian algebra of 
nonlocal symmetries. 

That option is not so apparent for the present system. However, we may make use of an 
extension of the work of Nucci et al [II] to obtain indirectly a representation of the complete 
symmetry group for the system . They, by a complex sequence of analyses and transformations 
facilitated by some recent theorems of Andriopoulos et al [I], were able to represent the complete 
symmetry group of a single Ermakov-Pinney equation. In the notation adopted here, (|9aj) say, 
is completely described by the three symmetries 

Ai = a t , A 2 = A 3 = \( [(l + ^)dt)d t . (16) 



2 

It is a simple matter to infer that the system © is completely specified by plus 



^=ilf)dt, A 5 = i(7(l + -l^)dt)^ (17) 



to give the five symmetries expected by the general theory [2] . 

We verify that Ai — A5 are a valid representation of the complete symmetry group by con- 

[21 [21 

sidering the actions of A^ — A 5 on the general two-dimensional system 

x = f(t,x,y,x,y) , 
V = 9(t,x,y,x,y) . 

The action of Ai removes t from both of / and g. The second extensions of A2 and A4 are 

= (/ ^2 dt ) d * ~ \ 9 * and A l 21 = (/ ^2 dt ) 9 t ~ \ d V- 

They remove x and y from / and g. The relevant parts of the second extensions of A3 and A5 
are 

A 3 2] = {-/ + -L} 5 * and *$ = {-g + l s }a B 
from which the recovery of © is evident. 

in itself is interesting as the three integrals Ji, I2 and ^3, respectively Ji, J2 and J3, are generalisations of the 
quadratic integrals of the one-dimensional free particle which for x = are |i 2 , x(tx — x) and \{tx — x) 2 . These 
integrals are obviously products of the two linear integrals x and tx — x. The nonlinear Ermakov system does 
not possess linear integrals. However, the quadratic integrals could be considered as the inner products of the 
nonconserved vectors i x = id.i + -&2, i y = ye.\ + -£2, j x = (tx-x)ai + -0,2 and j = (ty — y)e\ + -e.2, where 6i 

x y x y y 

and ej, i = 1, 2, are unit vectors. In further contrast to the free particle the inner products of the x vectors and 
the y vectors are not conserved. 
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4 The method of reduction of order 



The method of reduction of order was used very effectively by Nucci |38| to obtain the complete 
symmetry group for the Kepler Problem, the representation of which is mainly in terms of 
nonlocal symmetries, in terms of Lie point symmetries of the reduced system. The reader is 
referred to Nucci et al |39j for a detailed description of the method. Essentially the differential 
equations are written as a system of first-order ordinary differential equations and the order of 
the system is reduced by the removal of ignorable variables. Eventually one (or more) variable(s) 
is (are) eliminated to give at least one equation of the second order so that Lie's algorithm for the 
determination of Lie point symmetries can be implemented as a finite procedure. The symmetries 
obtained are then translated back to symmetries of the original system. The identification of 
the complete symmetry group is facilitated by the theorems on transitivity recently obtained by 
Andriopoulos et al [I]. 

The Ermakov System as presented in Q is not satisfactory for the application of the method 
of reduction of order since there is no preferred new variable after the time, as an ignorable 
coordinate, is eliminated. However, conversion of the system to a representation in terms of 
plane polar coordinates through x = r cos 9 and y = r sin 9 and recombination into radial and 
tangental components, videlicet 

r -r9 2 -r~ 3 (tan 9 + cot 0) 2 = 0, (18a) 
rO + 2r9 + ±r -3 (tan 9 + cot 9)' = 0, (18b) 



where the prime denotes differentiation with respect to 9, does provide a good basis for the 
reduction since ultimately one can envisage a polar representation of the orbit. 

We define the variables w\ = r, w 2 = 9, w% = r and ioj = 9 so that the system (fT8|) can be 
written in terms of the four first-order ordinary differential equations 

wi = w 3 , (19a) 

W2 = w 4 , (19b) 
2 a{w2) 

u>3 = wiw 4 H 3—, (19c) 



w 3 w A a'{w 2 ) 

m = -2 —4-, (19d) 

wi 2wt 



where we write (tan# + cot 9) 2 = a{9) to keep a certain compactness of notation. Obviously the 
system is autonomous, i.e. t is an ignorable variable. 

An obvious candidate as a new independent variable is w 2 which we write as y. The ratios 
of (|19a|) . ()19cj) and ()19d|) with ()19bl) give the reduced third-order system 

dwi w 3 

— — = — , (20a) 
ay W4 

— — = W1W4 H ^ — , (20b) 

ay w^W/i 

^± = _ 2 ^_^L. (20c) 

ay wi 2wjW4 

The next step is to introduce a second-order differential equation into the system so that Lie's 
algorithm becomes finite. The obvious variable for elimination is ti^ for from (|20aj) we have the 
simple relationship 103 = W4w[, where now the prime refers to y. 
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The system of second-order differential equation and first-order differential equation is 

W4W1 + w A w l = W1W4 H o — , (21a) 

wfw4 



, n W4W 1 ayy . , 

W i = -2— TT-i— • ( 21b ) 



) uhw[ _ a'(y) 
w\ 2w\w4 

At this stage one can apply Lie's algorithm, but some aesthetic reshaping of system (|2~Tj) is 
informative. Equation 1)2 lb|) can be written in the explicitly conservative form 

v! 2 = 0, (22) 

where 

u 2 := w\w\ + a(y). 

We recognise 1*2(2/) as related to the Ermakov Invariant appropriate to system @, but we must 
recall, as the explicit presence of the independent variable reminds us, that this is in terms of the 
independent variable y and not the original independent variable t. The distinction is technical 
since du 2 /dt = u' 2 y. 

We introduce the new dependent variable into (|21a|) and after some manipulation arrive at 

(1*2 — a) u'i — \a'u' x + 1*21*1 = (23) 

where we have introduced the new dependent variable u\ = \/w\. 

We have reduced the original fourth-order system to a third-order system comprising one 
second-order equation (|2"3*)) and one first-order equation (|22j) . The first-order equation has the 
nature of a conservation law. We recall that the same type of reduction has been demonstrated 
for the Kepler Problem, problems based upon it and the MIC-Z system. The question which 
naturally arises is one of the equivalence of these somewhat disparate systems. The resolution 
of the question is found in consideration of ()23|) . As a linear equation this can be reduced to 
that of the simple harmonic oscillator - the form of the second-order equation taken in the case 
of the reduction of the Kepler Problem - by means of a Kummer-Liouville transformation. The 
transformation of ()23|) to a nonautonomous oscillator requires simply a change in the dependent 
variable. However, to go from the nonautonomous oscillator to the simple harmonic oscillator 
requires a redefinition of the independent variable. In the usual area of application, i.e. the 
treatment of the time-dependent oscillator, the introduction of 'new time' is simply a rescaling 
of a variable without bound. Here the independent variable has a geometrical meaning and a 
rescaling of the polar angle changes the geometry. Consequently we must conclude that there 
is an essential difference between the two classes of problem despite the formal equivalence 
obtained under the procedure of the method of reduction of order. 

5 Symmetries of the original system 
and its complete symmetry group 

The Lie point symmetries of the system (|23|) . (1221) are given in composite form as 

r = (a + Pui)d y + {au\ + 7m + 6)d ui + ed U2 , (24) 

where 

* = f+2(^) ft < 25 »> 

7 = 70 + I(* + ^), (25b) 
2 \dy u 2 - a ) 



10 



P.G.L. Leach, A. Karasu (Kalkanh), M.C. Nucci and K. Andriopoulos 



70 is a constant, e is an arbitrary function of 112 and a, (3 and 5 solutions of the equations 

2" 



da 
dy 



4u, 



d 



u 2 — a dy \u 2 — a 



u 2 - a 



da 
dy 



ld_ 

+ 2dy 



4U2 



d 
+ t; 



1 



U2 — a dy \U2 — a 
d 

+ 



2(u2 — a) \ dy \U2 — a J U2 — a J ' 



1 

~ 4 
4a 



a 



u 2 - a 



a 



dp 1 a' 9/3 1 [ 2u 2 
<9y 2u 2 — a dy 2 \_u 2 — a 
d5 1 a' 55 U2 



U2 — a 



(3 = 0, 



+ 



<9y 2u2 — ady U2 — a 



0. 



(26a) 
(26b) 
(26c) 



We note that the equation for j3 is the adjoint of the equation for 5. 

Two of the relevant symmetries are the two solution symmetries of ()23|) given by the 5 
equation. We write these as 



sid Ul . 



S2d Ul , 



(27) 



where s± and S2 are two linearly independent solutions of the 5 equation. 

We firstly examine the effects of these two symmetries. The general pair of second-order 
equation and first-order equation is 

u'{ = f(y,ui,u 2 ,u' 1 ), 
u'2 = 9(y,ui,U2,u[). 

The application of gives 



df_ 

du\ 
dg 
du\ 



, df 

+ Sl d7 1 > 



= SlTT^ + Si 



1 dg 
du[ 



from the solutions of which the two equations in ()27|) now become 



- + F(y,s 1 u' l - s' 1 u' 1 ,u 2 ), 



si 



u' 2 = G(y,sini - s 1 u' 1 ,u 2 ), 

where F and G are arbitrary functions of the indicated arguments. 
After the application of T, 2 to ()28|) we have 

1 sis 2 - s'{s2 



(28) 



it--) 



— Ul + 
Si Si SlS 2 — S\S2 



siu^ - s[ui) + T{y,u 2 ), 



Q(y,U2) 



(29) 



Since 
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the system (|29|) becomes 



(u 2 - a)u[ - \au\ + u 2 u\ = F(y, u 2 ), 
u 2 = G(y,u 2 ). 

The two remaining symmetries required to specify the system come from the particular 
solution of (|26a|) and its consequence in ()25b|) . The particular solution of (|26a[) follows from 
that of the original equation (|23j) and the use of Laplace's method of variation of parameters. 
The symmetries can be written as 

£ 3 = A d y + C Q uid Ul + d U2 , 

£4 = A u 2 d y + (C u 2 + 7o)^i«9 ni + u 2 d U2 = u 2 £ 3 + 7 uA 2 , 



where 70 is a constant, Co is given by ()25b|) and 

-,'2 



A 







1 

w 



°\ \ a 2 



+ 



a 



+ 



4a 



0\0 2 / 0\0 2 



(u 2 - a) 2 (u 2 - a) 3 (u 2 - a) 2 

+ 7 no" 4" 



a 



(n 2 - a) 2 (n 2 - a) 3 (n 2 - a) 2 

a' 2 4a 
+ 7 ^ + 



dy 



(u 2 - a) 2 (u 2 - a) 3 (u 2 - a) 2 



dy 



In fact the complexity of the expression for the symmetry is not of importance. Rather it is the 
structure found in £3 and £4. 

A symmetry of the original system @ has the general form 

r = rd t + id x + r)d y 

in which the variable dependence in the coefficient functions r, £ and rj is not specified since 
for the symmetries of interest to us this is determined through the connection to the point 
symmetries listed in (|24j) . We follow the line of transformations through conversion to polar, 
i.e. x = rcos#, y = rsin# and 

T — ► rdt + (£ cos 8 + 77 sin 9) d r + - (— £ sin + 77 cos 0) <% 

r 

which for compactness of notation we write as 

r — > rd t + i?a r + ea , 

expression as the system of four first-order differential equations, i.e. r = w\, 6 = w 2 , r = W3, 
6 = W4 and 



rd t + R8 W1 + Od W2 + (R- rf)d m + (G - 0r)3, 



U'4 



the reduction of order which means the elimination of t and the use of w 2 = y as the independent 
variable giving the operator 

r — > Rd m + Qd y + (r- rf)d W3 + (e - ^f)a W4 . 

Finally we have the elimination of 77J3 and the introduction of the two new dependent variables 
u\ = l/wi and u 2 = w\w\ — a(y) so that 

r — ► (6 - (6 - 9f)a')d y - u\Rd Ul + (4wlwlR + 2wjw 4 (@ - Of))d U2 . 
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We are now in a position to express the symmetries (j27fl in terms of the original variables of 
our toy Ermakov system. After a modest amount of manipulation we obtain the symmetries 



Ai 



A, 



A, 




s 1 (9)rdt)d t + s 1 (e) (xd x + yd y ), 



s 2 (9)rdt)d t + s 2 {6) (xd x + yd y ), 



1 



9 2r 4 6» 2 
C y + x [A + 



2C 



dt\d t 



C x + y ( A + — -2 + 2C a 

2 (xy - xy) 



d x 



2 (xy - xy)' 

7 



+ 2C a 



A 4 = ((xy - xy) + a A 3 + — x ( xd x + yd y ) , 
\ J x z + y z 

A 5 = d t , 



(30) 



where the fifth symmetry is introduced as the symmetry responsible for the reduction of order by 
the change of independent variable from t to 6 and s\ and s 2 are linearly independent solutions 
of lf25]). 

In (|3Uj) we have a representation of the complete symmetry group of our original system jUJ). 
The algebra of the symmetries may be written as A\ ® {A 2 © s 2Ai} in the notation of the 
Mubarakzyanov classification scheme [3111341 1531 136 j . The first Abelian subalgebra comprises A5 
and the second abelian subalgebra the solution symmetries Aiand A2. 

Note that the four symmetries coming from the system (|22|). ()23|) are nonlocal only in the 
coefficient of dt- The symmetries Ai and A2 are point in the spatial variables as is the case with 
the nonlocal symmetries used by Krause to specify the Kepler Problem. The symmetries A3 
and A4 are generalised in the spatial variables. Given that the original system is autonomous, 
the concentration of the nonlocality in the coefficient of dt is not surprising since the integral 
does not enter of itself into the action of the symmetry on the system. 



6 Conclusion 

In this paper we have considered a very simple example of the class of problems known as 
Ermakov systems. The simplicity of the choice was deliberate for we were then able to concen- 
trate upon the important concepts rather than becoming enmired in distracting manipulations. 
Even still the conversion of the symmetries from their simple representations for the system (|22|) 
and ((2*3^1 to those given in (|3T)|) was not a completely trivial exercise. We have seen that a rep- 
resentation of the complete symmetry group of this toy Ermakov system cannot be achieved by 
means of Lie point symmetries. It is necessary to look to nonlocal symmetries. Nevertheless 
our toy system belongs to the class of superintegrable problems. It is also of some interest that 
by means of the method of reduction of order we were able to reduce the Ermakov system, 
fourth-order and nonlinear, to a third-order system of linear equations. This same reduction 
has been observed for other problems, but there seem to be some essential differences due to the 
geometric implications of the transformation of the independent variable. 

As a final remark we recall that the coupled system Q is superintegrable. Other superin- 
tegrable systems include the simple harmonic oscillator and the Kepler Problem. The former 
enjoys a richness of symmetry which is not lessened by the adjoining of a further dimension 
or two. The latter never has a sufficient number of Lie point symmetries for its complete 
specification be the dimension two or twenty-two. Yet it is the classic example of a natural su- 
perintegrable system. The Ermakov system considered here occupies a middle position in that 
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it has sufficient Lie point symmetries as a scalar system to be more than superintegrable and yet 
does not have the requisite number of point symmetries to be completely specifiable. Here we 
have detailed a strategy to reconcile the superintegrablity of the system and its ease of complete 
specification. The strategies for the resolution of particular problems may vary, but they have 
the same underlying basis which is to seek that representation in which the superintegrability 
is supported by an equally generous supply of Lie point symmetries for they are the class of 
symmetry with which to work is the easiest. 

The conventional representations of systems such as the cojoined Ermakov system considered 
here are not necessarily those for which the representation of the Lie algebraic properties are 
the simplest. 
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